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ABSTRACT 

We have used the very large JVAS/CLASS 8.4-GHz surveys of flat-spectrum radio 
sources to test the hypothesis that there is a systematic ahgnment of polarization 
position angle vectors on cosmological scales of the type claimed by Hutsemekers et 
al. (2005). The polarization position angles of 4290 sources with polarized flux density 
>1 mjy have been examined. They do not reveal large-scale alignments either as a 
whole or when split in half into high-redshift (z > 1.24) and low-redshift sub-samples. 
Nor do the radio sources which lie in the specific areas covered by Hutsemekers et al. 
(2005) show any significant effect. We have also looked at the position angles of parsec- 
scale jets derived from VLBI observations and again find no evidence for systematic 
alignments. Finally, we have investigated the correlation between the polarization 
position angle and those of the parsec-scale jets. As expected, we find that there is a 
tendency for the polarization angles to be perpendicular to the jet angles. However, 
the difference in jet and polarization position angles does not show any systematic 
trend in different parts of the sky. 
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1 INTRODUCTION 

In an isotropic universe the sizes of the largest coherent 
structures are expected to be limited by the time taken 
for the primordial density fluctuations to undergo gravita- 
tional collapse. In practice both observations of large-scale 
structures and the latest N-body simulations agree that the 
largest structures have scales of <100 Mpc. It is therefore ob- 
vious that any observation of coherent behaviour of objects 
on much larger scales than ~100 Mpc would have profound 
implications for our understanding of the nature of the Uni- 
verse. Over the years there have in fact been several reports 
of the detection of large-scale coherence, mostly based on 
measurements of the linear polarizations of active galaxies, 
and more recently based on observations of the cosmic mi- 
crowave background (CMB). 

Birch (1982) claimed to have found evidence for a ro- 
tating universe by comparing the radio polarization posi- 
tion angles for extended radio sources with the direction of 
elongation of radio sources. The claim was that there were 
systematic offsets between the two position angles and that 
these were systematically different in widely different part 
of the sky. Though attacked on statistical grounds at the 
time (e.g. Phinney & Webster 1983) it is not clear that this 



"Birch effect" has totally gone away (Kendall & Young 1984; 
Jain et al., 2003). An independent effect claimed by Nodland 
and Ralston (1997) was that there is a systematic rotation of 
plane of polarization over cosmological distances. But War- 
die et al (1997) refute their claim by pointing out that in 
individual radio sources the polarization position angles are 
strongly correlated with extended radio morphology, even 
at high redshift which excludes large cosmological rotations. 
Also the Nodland and Ralston claim has been attacked on 
the basis of statistics (e.g. Carroll and Field (1997), Eisen- 
stien and Bunn (1997), Loredo et al (1997)). 

More recently, Hutsemekers (1998), Hutsemekers & 
Lamy (2001) and Hutsemekers et al. (2005) have been ac- 
cumulating evidence that the linear polarizations of quasars 
in optical wavelengths are non-uniformly distributed, being 
systematically different near the North and South Galactic 
poles. They also find that the clustering of position angles 
is more evident if they divide their sample of 355 sources by 
redshift. They use this fact to argue against the effect being 
a result of the polarization of the quasar light as it passes 
through Galactic interstellar dust since this would affect all 
redshifts indiscriminately. 

Stimulated by these hints of interesting polarization ef- 
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fects, and reports of a preferred axis in CMB fluctuations 
(Schwarz et al., 2004; Land & Magueijo, 2005; Raeth et al., 
2007), we decided to examine some linear polarization mea- 
surements we had made with the Very Large Array (VLA) 
between 1990 and 1999 of a large sample of compact flat- 
spectrum radio sources. The JVAS (Jodrell Bank- VLA As- 
trometric Survey) programme from 1990-1992 (Patnaik et 
al. 1992, Browne et al. 1998, Wilkinson et al. 1998) was un- 
dertaken with the two main objectives of measuring accurate 
positions for sources to find interferometer phase calibrators 
and to identify gravitationally lensed sources. However, we 
also obtained data for the linear polarizations of the 2308 
sources. In addition, the later CLASS survey (1994-1999; 
Myers et al. 2003, Browne et al. 2003) targeted fainter com- 
pact radio sources and the majority of these observations 
have suitable calibrators available to derive degrees of polar- 
ization. The total size of the JVAS/CLASS surveys is 16503 
sources, and as we discuss in Jackson et al., (2007) (here- 
after Paper I), 4290 have polarized flux density >1 mjy. 
These 4290 sources we make use of for studies of possible 
polarization alignments. 

On its way to the observer polarized radiation passes 
both through the host galaxy of the AGN and the Galaxy 
and suffers Faraday rotation. Since it is dependent on the 
square of the wavelength, Faraday rotation is therefore a 
much more significant effect at radio than at optical wave- 
lengths. However, as we have shown in Paper I, total rotation 
measures are generally much too small to destroy informa- 
tion about the intrinsic position angles of polarization when 
the measurements are made at 8.4 GHz. 



2 INITIAL RESULTS AND REANALYSIS 

Our initial examination of JVAS for evidence of aligned radio 

polarizations made use of the original calibrations of the 
data done primarily with astrometry in mind and not the 
accurate measurement of source polarizations. Using these 
data we did initially find apparent systematic alignments so 
we decided on a careful re-calibration of both the JVAS and 
CLASS data, concentrating on polarizations. The results of 
this recalibration and analysis are presented in Paper L 

2.1 The data and analysis 

We make use of the 4290 sources presented in Paper I and 
full details of the calibration and analysis are given therein. 
See Fig. 1 for the Aitoff projection of the data. Here we only 
briefly address calibration and analysis issues which may be 
directly relevant to any search for systematic alignments of 
polarization position angles within the sample. These issues 
are: 

• Polarization residual calibration. Incorrect removal of 
instrumental polarization residuals can lead to a bias in the 
measured polarization angles. The JVAS/CLASS observa- 
tions were concentrated in different regions of sky in differ- 
ent observing runs and thus it is possible that poor, and 
epoch-dependent, instrumental residual removal could lead 
to apparent large scale alignments in different regions of sky. 
On the other hand, incorrect residual removal is very un- 
likely to mask real regions of aligned position angles as the 



errors would have to conspire in such a way as to randomize 
the observed position angles in regions of true alignment. 

The analysis presented in Paper I leads us to believe that 
systematic errors in the residual calibrations are at a level 
~ 0.3%. 

• Polarization position angle calibration. Errors in 
position-angle calibration would lead to position angles be- 
ing systematically wrong, perhaps by different amounts in 
different observing runs. We believe that errors in the po- 
sition angle calibration are at the level of < 10°. Position- 
angle calibration errors cannot produce areas of systematic 
alignments where none exist. They could, however, reduce 
the prominence of areas of real alignment if the errors were 
large enough and changed systematically across a region of 
real alignment. But at < 10° these errors are small and, in 
fact, comparable to the random errors on the majority of 
sources. Thus they are not a good reason for missing real 
alignments. 

• Data analysis - CLEAN bias. Even with perfect cali- 
bration it is possible for the analysis technique adopted to 
lead to systematic biases in the distribution of polarization 
position angles. In particular the CLEAN deconvolution pro- 
cess suffers from a tendency to bias flux densities towards 
zero by an amount that depends on the number of CLEAN 
iterations. For total intensity data this is relatively benign, 
but for the analysis of polarization data where Stokes Q and 
U maps are separately cleaned it can cause significant prob- 
lems when Q and U are combined to give a position angle. 
This is because Q and U can be negative and because polar- 
ized flux densities are generally low. Hence biases towards 
zero are proportionately more significant and result in a dis- 
proportionate number of sources having apparent position 
angles around zero, ±45° and ±90°. The issue of CLEAN 
bias and other systematic effects which may affect polariza- 
tion measurements are discussed more extensively in Battye 
et al. (in prep). It is to avoid the effects of CLEAN bias that 
the results we analyse in this paper have all been obtained 
by model fitting in the visibility plane without the use of 
CLEAN. 

2.2 Westerbork Data at 5 GHz 

JVAS and CLASS were made using the VLA which has feeds 
that produce opposite hands of circular polarization that 
can then be cross-correlated to produce linear polarizations 
(Stokes Q and U). We thought it useful to check the ac- 
curacy of the VLA polarization results by comparing them 
with measurements obtained in a different manner. Also, if 
the measurements were done at a different frequency this 
would give some Faraday rotation information. (See Section 
7) We therefore observed 340 sources from JVAS/CLASS 
lists with Westerbork Synthesis Radio Array (WSRT) at 
5 GHz. This uses linear-polarization feeds to obtain the po- 
larization information; thus the methodology is completely 
different from that of the VLA. 

All the sources observed were within < 5 < 40°, and 
had polarized fiux density, P > 4 mJy at 8.4 GHz. Not all 
the measurements were successful and we obtained useful 
polarization data for 336 objects. Fig. 2 compares position 
angles (PA) from JVAS/CLASS and WSRT data; the dif- 
ference clearly peaks at 0°. 

Full details of the observations, their calibration and 
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Figure 1. AitoflF plot for all the objects in the sample we used, from JVAS/CLASS samples. The lack of uniformity arises from the 
exclusion of low Galactic latitude sources from CLASS and the fact that in some CLASS regions we were not able to calibrate the 
polarizations to our satisfaction. 



analysis will be given elsewhere (Joshi ct al, in prep). Briefly, 
CTD 93 and 3C 286 were observed as flux and polariza- 
tion calibrators. The data analysis was performed using the 
NEWSTAR package. The target sources were assumed to 
be point-like and the values for Stokes Q, U (and V) were 
obtained by model fitting. 
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3 STATISTICAL ANALYSIS 

Our polarization data as described in Paper I have been 
analyzed for signs of statistical non-uniformity of the type 
claimed by Hutscmckcrs ct al (2005). In Paper I it was 
shown that the distribution of position angles for the sample 
treated as a whole was consistent with uniformity. However, 
this does not preclude the possibility of sub-areas in which 
there are locally strong alignments. Four tests have been 
performed. In the first two the sky has been divided into 
"tiles" defined using the HEALPix software package (Gorski 
et al., 2005) and the distribution within each tile has been 
tested for non-uniformity using a test on the histograms 
and by adding the polarization vectors to form a random 
walk. In the third test a search has been made for clusters 
of aligned position angles using a nearest neighbour analy- 
sis. In the last test we implemented one of the approaches 
used by Hutsemekers et al. All these tests applied to the 
4290 objects from JVAS and CLASS catalogues for which 
P > ImJy. 

3.1 Histogram method 

We have investigated a simple adaptation of the histogram 
test performed on whole dataset in Paper I which is designed 
to search for localized non-uniform behaviour. The data were 
split up into pixels defined by the HEALPix software and 
the of each of these histograms of position angles was 
computed on the assumption that the expected distribution 



Table 1. Results of the histogram test. 0pix is the approximate 
pixel size, riobj is the number of pixels containing an object. 
N{95%) and iV(99%) are the number of pixels which fail the test 
at 95% and 99% respectively, that is, those pixels for which the 

is greater than that found in 500 and 100 of the random 10000 re- 
alizations, respectively. P(95%) and P(99%) quantify the number 
of failures as a percentage of the pixels containing a source. 

was uniform. This value is then compared to the x^ found 
for 10000 realizations of sources at the same positions, but 
with PAs which are drawn from a uniform distribution. For 
the real data the number of pixels that have a x^ indicating 
that the distribution has failed the test for uniformity at 
different degrees of significance is noted. This is compared 
to the number one might expect for that number of pixels 

The size of each pixel is defined by the HEALPix vari- 
able nsidc, which gives the number of pixels in the whole 
sky ripix — 12ngid^, and hence the area of a pixel is f2pix = 
^pix = 47r/npix. We also need to choose a bin size for the 
histogram which is defined by nbin, the number of bins cov- 
ering the range —90° to 90°, which is varied in the analysis. 
Moreover, we have performed this test using two different co- 
ordinate systems (celestial and galactic) so as to investigate 
the effects of choice of pixel shape defined by HEALPix. 

The results of this test are presented in Table 1 for 
?^bin = 4 and the celestial coordinate system. For large pixel 
sizes (Spix ~ 58.5° and 29.3°) there is a slight excess of 
pixels which fail the test - one would expect around 5% of 
pixels to fail the test at 95% confidence, for example - and 
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Figure 2. The histogram of the differences between position angles measured using the WSRT and the VLA in the JVAS/CLASS 
surveys. 



substantially less failing than one expects for ^pix 3.7°. 
The latter result can be understood by realizing that there 
are, on average, less than 5 sources per pixel in this case. The 
results for large pixel sizes can be "explained" as being due 
to small number statistics. For 0pix = 14.7° and 7.3°, things 
seem compatible with there being no excess alignment, over 
and above that which is expected. 

This test is however not completely robust since it re- 
lies on a choice of ribin and the coordinate system. We have 
checked that the results presented in Table 1 are not partic- 
ularly sensitive to these choices by using ribin = 6, and also 
a galactic based coordinate system. Hence, we conclude that 
the test is compatible with the null hypothesis that the PAs 
are chosen from a uniform distribution. 

3.2 Random walk test 

We have performed another test using the HEALPix-defined 
pixelation scheme: for each source in the pixel we have con- 
structed a unit vector in the direction of the polarization 
{pi = (sin 20, cos 20)) and these are then added together 
to form a random walk. Again the results are compared to 
10000 realizations whose PAs are drawn from a uniform dis- 
tribution. In contrast to the histogram test, this docs not 
require the choice of ribin; it also seems to work better when 
there are smaller numbers of sources in the pixels. It does, 
however, still depend on a choice of coordinate system. 

The results of this test are presented in Table 2. We see 
that again that there is a slight excess of pixels which fail 
the test for large pixel sizes, but that the number failing the 
test appears compatible with uniformity for lower values of 

3.3 Nearest Neighbour Method 

As a simple global test to sec if the PA arc aligned we used 
the following method: for every object, the n nearest objects 
are checked for alignment. If the object has PA within AO 
of the object under consideration, it is considered aligned. 
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Table 2. Results from the random walk test. The columns are 
the same as in Table 1 



The choice of the appropriate value of n depends on the size 
of the sample, since the probability of finding n alignments 
by chance should be low. It also depends upon the scale of 
any real clustering in the data which is, of course, unknown. 
However, if the regions over which alignments occur contain 
more than n objects, the test will pick out many objects in 
the same area all of which will have significantly more rear 
neighbours having aligned polarizations than expected by 
chance and real structures will not be missed. We initially 
choose a value of n = 25 but have explored other values 
up to a maximum of n = 250. We typically used AO = 45° 
since all the errors on the PA measurements are significantly 
smaller than this angle. 

For the whole sample a histogram showing the number 
of objects having N out of 25 aligned nearest neighbours is 
plotted in Fig 3. The error bars have been derived from 5000 
random realizations where, keeping positions of the objects 
the same, the sky was populated with sources having random 
PAs and the nearest neighbour test performed on each real- 
ization thus producing 5000 histograms like the one shown 
in Fig 3. For each bin in the histograms the distribution of 
values was obtained and fitted by a Gaussian. The plotted 
points and error bars are, respectively, the means and stan- 
dard deviations obtained from the Gaussian fits. It is clear 
that, given the error bars, the real histogram is a reasonably 
good fit to that expected for a sample with randomly orien- 
tated position angles; the for the real distribution when 
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Figure 3. Compajrison of nearest neighbour distribution for the 4290 JVAS/CLASS objects with that expected for a randomly oriented 
sample obtained from 5000 random realizations. The solid lined-histogram is for the real data and dashed data with error bars is from 
5000 random samples, this distribution is 20.1. The right hand panel shows the distribution of obtained by comparing each 

random realization with the averaged of all the random ones. The line represents x^ for the observed sky. 



compared with the average random distribution is 20.1. The 
appropriate number of degrees of freedom is not obvious 
because the numbers in the histogram bins are not statisti- 
cally independent. In Fig. 3b wo show the distribution of x^, 
obtained from the random realizations, by comparing each 
with the mean of the 5000, with the line marking the 
for the real sky. The probability of finding a equal to or 
greater than the observed one is 45%. 



3.4 Hutsemekers' test 

Hutsemekers (1998) used a dedicated statistical test to iden- 
tify if the polarization vectors arc aligned in 2D or 3D space. 
It is more sophisticated than our own nearest neighbour test 
in that it gives extra weight to aligned objects according 
to their degree of alignment and distance apart. (See Hut- 
semekers (1998), Section 5.1 for the complete description.) 
As in our irnplcrricntation of our nearest neighbour tost, 
alignments were searched for within ~45° in groups of 25 
objects. We have implemented the Hutsemekers (1998) test 
on ,1VAS/CLASS data in 2D (since we have incomplete red- 
shift information). In this test, a parameter, S is calculated 
which depends upon the spatial distribution of objects, their 
PAs and the distribution of these. We find S = 68.6 for the 
4290 sources in the JVAS/CLASS sample. To test whether 
or not this value of S indicates statistically significant clus- 
tering we generated 5000 random samples. For each sample 
the positions of the objects are same as those observed but 
PAs axe randomized. The S distribution obtained from the 
random samples is shown in figure 4. Using this distribution 
wo find that the probability of finding S > 68.6 occurring by 
chance is 19%. 



4 DIVISION OF THE SAMPLE INTO 
REDSHIFT BINS 

Hutsemekers et al. have found that the statistical signifi- 
cance of their alignments axe enhanced if they divide their 
sample into low and high redshift bins. We do likewise for 



the JVAS/CLASS results in order to make as close as pos- 
sible comparison with results of Hutsemekers et al. Out of 
4290 objects, redshift information is available for 1273 ob- 
jects. The redshifts arc between 0.03 and 4.72, with median 
redshift of 1.24. Again, the same nearest neighbour test was 
performed on the two samples, and on 5000 random realizar 
tions keeping the object positions the same, and assigning 
the position angles randomly from a uniform distribution. 
The histograms with error bars arc shown in Fig. 5, together 
with the distribution of x^- For the real sky, x^ values are 
23.8 and 18.8 for low and high redshift regions, respectively. 
Values of these or greater have probabilities of occxirrcncc 
as almost 1 in 2 and 1 in 4, in low and high redshift regions, 
respectively. 



5 AN EXAMINATION OF THE 

HUTSEMEKERS ET AL. REGIONS 

Hutsemekers et al. claim that in different regions of sky 
quasar optical polarization position vectors are strongly 
aligned. In Fig 6 we plot the distribution of radio polar- 
ization position angles of our sources that are found within 
the regions identified by Hutsemekers et al. They defined 
regions Al and A3 as, 168° < a < 218° and S < 50°, and 
320° < a < 360° and 5 < 50° , respectively. We do not have 
any radio data for 5 < 0° so we do not have exactly the 
same area coverage as they do. No obvious alignments are 
visible on the sky or in the histograms of position angles also 
shown in Fig. 6. We have also performed our nearest neigh- 
bour test on the observed data in those regions and again 
produced 5000 random realizations as described above. The 
resulting histograms are shown in Fig. 7. The two regions 
have x^ of 34.1 and 28.3 respectively, with the probability 
of occurrence of a x^ greater than or equal to these values 
being almost 1 in 2 and 1 in 5, respectively. In contrast, our 
test applied to the optical data, showed a high significance 
detection of non-uniformity. The probabilities of the optical 
position angles being uniformly distributed are < 0.7% and 
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Figure 4. The distribution of the parameter S for 5000 random samples. For the observed data, S is 68.6 and is marked by a line in the 
figure. 
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Figure 5. Histograms for low and high rcdshift samples. The panel on the top left shows histogram for objects with low rcdshift, the 
one at the bottom left is for high redshift objects. The solid line represents observed the data and dashed line is for random objects. The 
reduced x ^ for these two histograms are 0.26 and 0.29 respectively. 



12%, for regions Al and A3, are respectively^ (See Fig. 7, 
panels at the bottom for region Al.) 

Finally we compare the radio and optical position an- 



It should be noted that Hutsemkekers et al. preferentially tar- 
geted quasars in redshift slices where alignments had been pre- 
viously suspected and thus the non-uniformity might be more 
prominent than in the radio where no redshift targeting was done. 



gles for objects common to JVAS/CLASS and Hutsemekers 
et al. The areas do not overlap completely as Hutsemekers 
et al., have analyzed the data for alignments near northern 
and southern galactic poles, whereEis JVAS/CLASS surveys 
covered all RA ranges but only declination > 0° . We find 52 
of their sources for which we have reliable polarization mea- 
surements at 8.4 GHz. (See appendix for the data.) Figure 8 
shows the distribution of PA differences in optical and radio. 
There is no obvious correlation between the optical and ra- 
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Figure 6. The left top panel shows flat-sky projection of the 8.4 GHz PAs in Hutsemekers et al. region Al, and the right hand panel 
shows histogram of the PAs. The bottom panels are the same for Hutsemekers et al. region A3. 



dio polarization position angles. This is perhaps not surpris- 
ing since correlations between radio and optical polarization 
position angles have been searched for in the past but none 
conclusively found (Lister & Smith (2000)). In summary, we 
find no evidence to suggest that there is anything special 
about the radio polarization properties of the objects found 
in the Hutsemekers et al. regions. 



6 COMPARISON BETWEEN POLARIZATION 
POSITION ANGLES AND THE AXES OF 
PARSEC-SCALE JETS 

Birch (1982) claimed that there were systematic offsets be- 
tween radio polarization position angles and the direction of 
elongation of radio sources and that these were systemati- 
cally correlated over the large regions. We have undertaken 
a similar investigation using our polarization position angle 
measurements and jet axes derived from VLBI maps. We 
use information on sub-kiloparsec scales rather than much 
larger scale information used by Birch because the majority 
of JVAS/CLASS objects, as a result of the spectral index se- 
lection used to define the sample, are compact with little ob- 
servational information available on their structures on the 
scales of 10s to 100s of kiloparsecs. However, there is plenty 
of parsec-scale information. Using data from the literature 



(Beasley et al., 2002; Fomalont et al., 2003; Henstock et al., 
1995; Kovalev et al., 2005; Petrov et al., 2005; Petrov et al., 
2005; Kovalev et al., 2007; Polatidis et al., 1995; Taylor et 
al., 1994; Thakkar et al., 1995; Xu et al., 1995; Zensus et 
al., 2002) and that available in web-based archives we have 
estimated jet position angles for 1565 sources. There are 842 
of these for which we also have polarization position angles. 
(See online data.) The jet position angle is defined to be in 
the range —180° to 180° measured from North through East. 
The position angles were measured by eye aud we estimate 
the likely uncertainty of these measurements to be ±10°. 

The first thing to note is that the distribution of jet po- 
sition angles is not uniform (see Fig. 9a). The two peaks at ± 
90° arise because typical VLBI arrays have more resolution 
in the East- West direction than they do in the North-South 
direction. There is, therefore, a higher probability of being 
able to detect a measurable jet in the East- West direction 
than in the North-South. We note, however, that the dis- 
tributions of jet angles should not depend on the areas of 
sky in which they are measured. The non-uniform overall 
distribution should not severely afi^ect local statistics such 
as those employed in the earlier analysis of polarization po- 
sition angles. 

In order to investigate local correlations, we have done 
two things: we have examined the distribution of position 
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Figure 7. The left top panel shows the histogram from nearest neighbour test for region Al, with error bars from 5000 random 
realizations. The right top panel shows the distribution for random realizations with the line showing that for the real sky in region 
Al. The middle panels show the same, for region A3, for the region Al is 34.1 and 28.3 for the region A3. The probability of finding 
this or more is 52% for region Al and that for region A3 18% respectively. The panels at the bottom are for region Al again, but 
show the results for the optical data. For these x^ is 245 with probability of 0.7% under the hypthesis that the distribution of position 
angles is uniform. 



angles to see if there are any systematic deviations from 
uniformity in different parts of the sky and wo have looked 
at the polarization position angle/VLBI jet angle difference 
for any similar effects. 

6.1 The distribution of jet position angles 

We have performed the histogram and random walk tests 
on the data using nbins = 4 and the celestial coordinate 
system. The results are presented in Table 3. It is clear that 



for 6pix < 29.3° the results are compatible with uniformity 
and the case of ^pix = 58.6° is dominated by small number 
statistics. 

We also have performed the nearest neighbour test on 
the jet angle data to look for regions of sky where there 
might be clusters of sources with aligned position angles. No 
obvious regions were found. From these three tests, there- 
fore, we conclude that the jet angles show no tendency to 
align in different regions of sky. 



7 
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Figure 8. The PA difference distribution for 52 common sources in JVAS/CLASS survey and Hutsemekers et al (2005) sample of quasars. 
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Figure 9. a) The left hand panel shows jet-angle distribution of 1565 sources. Two peciks at ± 90° stand out clearly, b) On the right 
hand is the Aitoff projection of distribution of jet PA on the sky. 



6.2 Correlations between jet and polarization 
position angles 



regions having large numbers of either negative or positive 
position angle differences. 



The histogram of the difference between jet PA and polar- 
ization PA is shown in Fig. 10. It is clear that there is a 
peak near 90° with a less significant peak around 0° . These 
90° and 0° peaks are consistent with previous results (e.g. 
Gabuzda ct al., 1994; Helmboldt, et al., 2007; Lister and 
Smith 2000; Pollack et al., 2002). Such a correlation is also 
expected on astrophysical grounds since the magnetic field 
vectors of the emitted radiation are known to align along 
the local jet direction or, less frequently, perpendicular to 
it (e.g. Gabuzda et al., 1992). The lack of uniformity in the 
polarization/jet position angle distribution, which survives 
the effects of resolution and selection involved in the con- 
struction of our sample, means that one needs to be careful 
when looking for signatures of systematic alignments. How- 
ever, an obvious test is to look at the sign of the position 
angle difference for any evidence of an effect similar to that 
found by Birch (1982). We have done this using the nearest 
neighbour test where we define a 'success' as a neighbour 
having the same sign as the source in question and a failure, 
one having the opposite sign. The results of the test are dis- 
played in histogram form in Fig 11. There is no evidence for 



7 DISCUSSION 

None of the statistical tests we have performed (see Section 
3) indicate any areas in which there are large-scale align- 
ments of PAs measured at 8.4GHz, in contrast to what has 
been found in the optical wavelengths by Hutsemekers et al. 
Could this be due to some error in the radio results aris- 
ing from instrumental bias in PA measurement, or a data 
analysis error? We think this is highly unlikely because the 
data have been carefully re-analyzed paying particular at- 
tention to the polarization calibration and the extraction 
of the polarization parameters. External comparisons with 
other data, including our own Westerbork observations, con- 
vince us that CLASS polarization angles are reliable (See 
Section 2.2.) 

We can think of two generic reasons why there might 
be alignments found in optical polarization position angles 
and not the radio ones: 

• The physical mechanism that aligns the optical polar- 
izations may be frequency dependent and not apply to the 
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Spix/deg 


"obj 


JVhi,t(95%) 




JVrand(95%) 


•Prand(95%) 


58.5 


12 


1 


8.3 


2 


16.0 


29.3 


44 





0.0 


1 


2.3 


14.7 


156 


4 


2.6 


6 


3.9 


7.3 


549 


2 


0.36 


17 


3.1 



Table 3. Results of using the histogram and random walk tests on the jet position angle data. The first two results are the number and 
percentage of pixels which fail the histogram test at 95% and the second two are for the random walk test. 
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Figure 10. The histogram of difference between VLBI jet PA and polarization PA. 



radio. This could either apply to the production of the ra- 
diation or any polarization imposed during propagation. 

• The optical results could be misleading because of ei- 
ther statistical fluke or the polarizations are not primarily 
extragalactic in origin. 

We first consider physical mechanisms. We have checked 
the optical sample of 355 quasars using Veron & Veron 
(2007) and the NASA Extragalactic Database (NED) and 
find that ~72% are radio loud and most have flat radio spec- 
tra. Having a flat radio spectrum suggests that they are 
blaaaxs and in blazars the optical and radio emission is pre- 
dominantly synchrotron in origin. Thus for the majority of 
the radio and optical objects the radiation mechanism is the 
same and this argues against something intrinsic producing 
alignments detectable at optical wavelengths and not at ra- 
dio wavelengths^. 

For radio-quiet quasars, the source of intrinsic optical 
polarization is likely to be scattering (Stockman, Moore & 
Angel (1984), Berriman et al (1990)). Perhaps significantly, 
the degree of alignment in radio-loud quasars and radio- 
quiet quasars cannot be distinguished statistically, though it 
should be noted that the radio-quiet sample is small. Thus, 
assuming the optical result is of cosmological significance, 
the more plausible explanation would be a propagation ef- 
fect that polarizes optical emission and not the radio, or 
perhaps one that destroys the radio alignments. An obvious 
mechanism to destroy radio alignments is Faraday rotation. 

^ Despite the common radiation mechanism we and other authors 
do not find any clear correlation between the radio and optical 
position angles. We have 52 objects in common with Hutsemekers 
et al (see Section 5) 



However, rotation measures would have to be of the order of 
thousands of radm~^ and this possibility can also be ruled 
out because: 

• The rotation measures of extragalactic sources in gen- 
eral are normally quite low, often less than a few tens of 
radm~^ (e.g. Rudrnck & Jones (1983), but see also Zavala 
& Taylor (2004) who measure very high rotation measure 
but on parsec scales). For JVAS/CLASS the random subset 
of objects for which we have measurements done at 8.4 GHz, 
5 GHz (WSRT), and 1.4 GHz (NVSS) the rotation measures 
arc^ around tens of radm^^ or less, rather than the order of 
thousands. (Joshi et al, in prep). 

• There is a correlation between the position angles of 
radio jots and the corresponding magnetic field directions 
inferred from radio polarization measurements (Wills et al., 
(1992); Visvanathan & Wills (1998); Fig. 10 of this paper). 
Such a correlation would not be seen if the radio polariza- 
tions were significantly randomized by Faraday rotation. 

We conclude that intrinsic mechanisms or Faraday ro- 
tation cannot account for the lack of radio alignments and 
therefore the most viable explanations are that either there 
is some, as yet uncertain, propagation effect that works on 
the optical photons and not the radio, or that the optical 
results have may not been interpreted correctly. This might 
be a result of the relatively small numbers alone or from the 
combination of small numbers and any residual biases in 
the data, perhaps caused by the effects of Galactic interstel- 
lar polarization. One suggestion discussed in Hutsemekers 
et al. (2005) is that radiation propagating over cosmolog- 
ical distances can be polarized by light weakly interacting 
pseudoscalar (or scalar) particles (Das et al, 2005) and that 



Aligned Radio Polarizations? 11 




No of successes out of 25 



Figure 11. The left hand panel shows the results of the nearest neighbour test on differences between jet angles and polarization position 
angles. The right hand panel shows distribution of with the real data having = 13.5, with probability of 69% of occurring under 
the hypothesis of a uniform distribution. 



this propagation effect can be frequency-dependent. How- 
ever, given the profound impUcations of a non-Galactic as- 
trophysical origin for the observed optical result, and the ab- 
sence of evidence for the equivalent effect seen in the radio, 
wo remain unconvinced by the conclusions of Hutsemekers 
et al. 

Finally, since there are models which predict that po- 
larization PA could suffer small (~ 1°) systematic rotations 
on large angular scales (e.g. Skrotskii rotation produced by 
vector perturbations of a flat Friedmann-Robcrtson- Walker 
background described in Morales & Saes, (2007)), we briefly 
address the question of whether or not in the future it might 
be possible to detect such small deviations from random 
alignments using radio polarization data. We suggest that 
the most promising approach would be to use the difference 
in PA between that of the polarization and that of parsec- 
scale jets measured with VLBI. The two are observed to 
be correlated with a dispersion of ^SO"'* and only the po- 
larization angles should be affected by the Skrotskii rota- 
tion. Approximately a thousand objects would be required 
to constrain the mean position angle at the degree level 
(30°/Vl000). To detect an effect, several regions containing 
the order of 1000 objects would be required. Thus, assum- 
ing no systematic errors, in samples of ~ 10* objects one 
would begin to be sensitive to alignments at the predicted 
levels. Already the number of objects with both polarization 
position angle measurements and jet measurements are sig- 
nificantly iu excess of a thousand, so such tests of this type 
are beginning to be fesisible and will be very powerful in the 
SKA era. 
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Table 4. The 52 common objects in JVAS/CLASS surveys and Hutsemekers et al (2005) sample. Position information is from 

JVAS/CLASS data. 1(8.4) stands for total flux at 8.4 GHz in mjy, p%(8.4) stands for percentage polarized flux at 8.4 GHz, p%(op)stands 
for percentage polarized flux in optical and PA(op) and PA(8.4) stand for polarization position angle in optical and at 8.4 GHz respec- 
tively. PA diff = PA(op)-PA(8.4). 



